Magneto-transport in Disordered Graphene: 
from Weak Localization to Strong Localization 



O 

(N 



in 



I 

O 

o 



> 

On 



m 
o 
o 



X 




J. Moser/ H. Tao/'B S. Roche/- ^ F. Alsina/ C. M. Sotomayor Torres,^-^ and A. Bachtold^ 

^CIN2 (CSIC-ICN) Barcelona, Campus UAB, E-08193 Bellaterra, Spam 
^CEA, IN AC, SP2M, L.sim, 17 avenue des Martyrs, 38054 Grenoble (France) 
■'^ICREA, 08010 Barcelona, Spain 

We present a magneto-transport study of graphene samples into which a mild disorder was in- 
troduced by exposure to ozone. Unlike the conductivity of pristine graphene, the conductivity of 
graphene samples exposed to ozone becomes very sensitive to temperature: it decreases by more 
than 3 orders of magnitude between 100 K and 1 K. By varying either an external gate voltage or 
temperature, we continuously tune the transport properties from the weak to the strong localization 
regime. We show that the transition occurs as the phase coherence length becomes comparable to 
the localization length. We also highlight the important role of disorder-enhanced electron-electron 
interaction on the resistivity. 



INTRODUCTION 

Quantum interference phenomena in graphene are of 
fundamental interest [ij. A case in point is the localiza- 
tion of charges, which is a manifestation of two important 
properties of this material: First, graphene hosts chiral 
Dirac fermions. Second, these fermions reside in two in- 
equivalent valleys at the K and K' points of the first 
Brillouin zone. Traveling paths that are relevant to lo- 
calization phenomena are phase coherent closed loops. 
Because of its chirality, a Dirac fermion residing in a 
given valley acquires a phase of tt upon completion of 
one loop, which gives rise to destructive interference with 
its time-reversed counterpart. Chirality therefore low- 
ers the probability for returning paths, and favors weak 
antilocalization. Restoring constructive interferences re- 
quires inter-valley scattering events (fermions in the K 
and K' valleys have opposite chiralities). This in turn 
favors weak localization. 

These quantum interference effects have been actively 
studied both theoretically [l-Q and experimentally (3- 
12|. It has been predicted that weak localization correc- 
tion to the semi-classical (Drude) conductivity will dom- 
inate the weak antilocalization correction as the temper- 
ature is lowered, driving graphene to the strong localiza- 
tion regime Surprisingly however, these corrections 
measured in graphene samples have remained modest, 
even at milliKelvin temperatures. 

A clear strategy to study the transition between weak 
localization and strong localization is to enhance interval- 
ley scattering. This can be achieved by introducing short- 
range scatterers, such as weak point disorder or lattice 

H-Eli. 



FIG. 1: (a) Atomic Force Microscopy image of the device 
contacted with 4 electrodes. The white contour highlights 
the sample shape. 



Approaches to drive the metallic phase of graphene 
to an insulator with an energy band gap have also been 



explored. In the case of graphane [17[, sp bonds were 
partially transformed into sp^ by hydrogenation. There, 
the resistivity was found to diverge at low temperatures, 
in accordance with the two-dimensional variable range 
hopping model. The measurements were interpreted as 
the result of a modified graphene that consists of two 
phases, regions with sp^ hybridization interspersed with 



gp2 regions 17|. Similar results were obtained with ox- 



defects that result in midgap states |13l - ll5| . Recently, 
defect scattering centers were introduced in graphene us- 
ing Ne and He ion irradiation, but the conductivity at the 
Dirac point remained above e^//i even down to cryogenic 
temperature {e^ /h is the conductivity value for which the 
weak localization regime is expected to cross over to the 
strong localization one [l6|). 



idized graphene 3 3|- Another work reported trans- 
port measurements deep in the strong localization regime 
using graphene modified with hydrogen atoms f20|. 

Overall, these transport studies on intentionally dis- 
ordered graphene focused either on the metallic regime 
[lij or on the deep localization regime, where the mate- 
rial behaves as an insulator [17H20| . The possible tran- 
sition between weak localization and strong localization 
remains to be explored. 

In this work, we report on a detailed magneto- 
transport study unveiling quantum localization effects 
in disordered graphene. Disorder is created by exposing 
graphene to ozone, which introduces sp^-type defects. By 
varying either an external gate voltage or temperature, 
we continuously tune the transport properties from the 
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FIG. 2: Raman spectrum of the device shown in Fig. 1 before 
(top) and after (bottom) exposure to ozone. G and 2D modes 
are shown; defects created by ozone reveal themselves as a 
strong D mode. Both panels have the same intensity axis. 



weak localization to the strong localization regime. We 
show that the transition occurs as the phase coherence 
length becomes comparable to the localization length. In 
addition, we evidence a marked contribution of electron- 
electron interaction to the resistivity. (This latter effect 
is a correction to the density of states and its origin dif- 
fers from the one of weak localization, which stems from 
a modification of the diffusion constant 



21|) 



EXPERIMENTAL RESULTS 

Before creating defects in graphene, we first fabri- 
cate high-quality devices using conventional nanofabri- 
cation techniques (Fig. 1) [22]. We mechanically exfoli- 
ate graphene from a flake of Kish graphite on a Si wafer 
coated with 300 nm of thermal oxide. We pattern Cr/Au 
electrodes in a four-point configuration using electron- 
beam lithography. We carry out Raman spectroscopy 
[23I and low temperature transport measurements to ver- 
ify that single layer graphene sheets are of good quality. 
For the device discussed in the paper, the D peak is ab- 
sent before ozone treatment (Fig. 2, upper panel). In 
addition, the mobility is 5500 cm^/Vs, and the conduc- 
tivity at the Dirac point reaches 5 e'^/h and is temper- 
ature independent down to liquid helium temperatures 
(not shown). 

We introduce defect using an ozone treatment, which 
is a chemically reactive process known to alter the under- 
lying sp^ network of graphitic systems 0. Specifically, 
we first clean graphene samples by placing them in a flow 
of Ar/Il2 gas at 300°C for 3 hours. We then expose the 
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FIG. 3: (a) 4-point resistance and (b) 4-point conductivity as 
a function of backgate bias Vg. Temperatures are T — 1.7, 3, 
5, 10, 20, 40 and 100 K. 



samples to ozone in a Novascan ozone chamber, where 
ozone is produced by ultra-violet irradiation of O2 gas 
(^ 7 min, 4 atm). The appearance of the D peak in the 
Raman spectrum (Fig. 2, lower panel) and the decrease 
of the mobility down to 390 cm^/Vs (see below) signal 
the creation of additional defects. However, we make 
sure that the ozone treatment is mild enough to preserve 
the crystalline integrity of graphene, as evidenced by the 
presence of the G peak and a well defined 2D peak in 
Raman spectroscopy. In addition, the elastic mean-free 
path of electrons estimated from the Drude conductivity 
is at least 3 nm (see below), which is more than one order 
of magnitude larger than the C-C bond length. 

This ozone treatment has a large impact on the trans- 
port properties of graphene. The conductivity a becomes 
very sensitive to temperature T (Fig. 3). Even though 
the conductivity at the Dirac point remains larger than 
e^/h at the highest temperature (100 K), it is reduced 
by more than 3 orders of magnitude at 1.7 K. The device 
behaves as an insulator, at least at low temperature and 
in the vicinity of the Dirac point (where the conductivity 
(T as a function of backgate bias Vg is lowest, see inset to 
Fig. 3). 

To reveal the contributions of quantum interferences, 
we explore the transport properties in the presence of a 
magnetic field B (Fig. 4a-f). In all cases, the magnetore- 
sistance is negative, and the resistivity changes with B 
by an amount that strongly depends on Vg. This change 
becomes increasingly large as the density approaches the 
Dirac point {Vg = 36 V, Fig. 4e) whereas it remains mod- 
erate at high charge density n {Vg = —40 V, Fig. 4a). 
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FIG. 4: Resistivity as a function of magnetic field, for T — 1.7, 5, 15, 50 and 100 K. Fits to eq. (1) are shown by dashed hnes 
away from the Dirac point and for conductivity values larger than /h. Inset to Fig. 4e displays the resistivity at the Dirac 
point at 5, 15, 50 and 100 K (top to bottom). 



DISCUSSION 

Both the low temperature insulating behavior [3, H] 
and the sharp D peak in Raman spectra show that the 
ozone treatment introduces significant intervalley scat- 
tering. A way to quantitatively compare intervalley 
and intravalley scattering rates is to examine the Vg 
dependence of the conductivity. Assuming only weak 
point disorder and charged impurity disorder, we ob- 
tain that the intervalley and the intravalley scattering 
times are comparable, about 10 fs for Vg = —40 V. 
Conversely, before ozone treatment the intervalley scat- 
tering time is 300 fs and the intravalley scattering time 
is 100 fs at the same n. To determine these scattering 
times, we consider weak point disorder for the interval- 



ley scattering time Tu 



(Tsr • /i/(2e^ • vpy/T^n) and 



charged-impurity disorder for the intravalley scattering 
time Tintra = h/{2e ■ vf^/tt) ■ [i^ - \jTi. The conductivity (Tsr 
[due to short range scattering] and the mobility [i can be 
extracted from the cr(V^) measurement taken at 100 K 
(assuming that localization effects are vanishingly small 



at higher temperature) using a = Ug 



(e/in)- 



We get CTgr = 5- 10~^ S and ^ = 5500 cm^/Vs for the pris- 
tine graphene, and Gsx — 4.1-10^'' S and [i = 390 cm^/Vs 
after ozone treatment. We note that lattice defects re- 
sulting in midgap states were recently identified as a new 



source of intervalley scattering 16| . The latter scattering 
results in a linear n dependence of a so its contribution to 
the conductivity cannot be discriminated from the one of 
charged-impurity disorder in a cr(V^) measurement. As 
such, it can modify rinter and Tintra and the intervalley 
scattering time of 10 fsec is an upper bound. 

The magnetoresistance measurements at high n can 
be well described by the weak localization (WL) theory 
developed for graphene @. The correction to the semi- 
classical (Drude) conductivity reads 
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with F{z) = In z + 'tp(0.5 + z^^) where ?/; is the digamma 
function, Tg^ = ieDB/h, and — L^/D the phase 
coherence time. To compare with the experiment, we 
let the diffusion constant D = 0.5vp ■ {r^lg,. + ^jntra)"^ 
and consider only weak point disorder for intervalley 
scattering and charged-impurity disorder for intravalley 
scattering. Accordingly, the phase coherence time 
is the only fitting parameter necessary. As illustrated 
in Fig. 4(a,b,c) (where a{B = 0) > e'^/h), we find a 
good agreement between experiments and theory. A 
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FIG. 5: (a) Temperature dependence of the resistance at the 
Dirac point, (b) At low temperature, the sample resistance is 
fit to the two-dimensional variable-range hopping model for 
several gate biases in the vicinity of the Dirac point. 



satisfactory agreement is also obtained by comparing 
measurements to WL predictions for conventional two- 
dimensional metals, which, moreover, yields the same 
phase coherence time. As for the magnetoresistance mea- 
surements at lower n (Fig. 4(d,e,f)), the resistivity can 
change with B hy a large amount. Comparing the mea- 
surements to theory is however difficult at this stage. 
More measurements at low temperature will be needed to 
discern among the various predicted dependencies p6| . 

We now turn our attention to the temperature depen- 
dence of the resistance. The resistance at the Dirac point 
diverges at low temperature, as expected for an insulat- 
ing regime (Fig. 5a). Fig. 5b shows that the temperature- 
dependent resistance is consistent with two-dimensional 
variable-range hopping, R cx exp((To/r)i/3). Even 
though the measurement could also be described by a 
simple thermal activation behavior, for the time being we 
will restrict our analysis to a variable range hopping sce- 
nario, which is the conventional mechanism to describe 
low t emp erature conduction in strongly disordered mate- 
rials |17h19| . |27| . This allows us to extract the localization 
length ^vRH from the fitting parameter Tq 



^VRH 



13.8 



(2) 



where p is the density of states of graphene. To esti- 
mate p, we assume that n — Cg^/V^ + Vq /e to take the 



fluctuations of the Dirac point with respect to the Fermi 
energy into account [11]. Here, Cg = 1.15 • 10"^ F/m^ is 
the backgate capacitance [2^. Letting Vq ~ 8 Y (which 
roughly accounts for the smoothing of the a vs Vg curve 
around the Dirac point at the highest temperature), we 
obtain £,vrh = 15 nm at the Dirac point. Away from the 
Dirac point, we find that ^vrh increases upon increas- 
ing n (Fig. 6a). As a comparison, we can evaluate the 
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FIG. 6: (a) Characteristic lengths as a function of backgate 
bias Vg. Z/e is the elastic mean free path. Coherence lengths 
at T = 1.7, 3, 5, 15, 25, 50, 70 and 100 K (top to bottom) 
are shown. (,vrh is the localization length estimated from the 
variable range hopping model (eq. 2) and is the localization 
length estimated from the scaling theory (eq. 3). (b) Phase 
diagram (temperature T, charge density n) showing the semi- 
metallic regime where weak localization (WL) is observed and 
the insulator regime where strong localization (SL) prevails. 
The top axis shows the corresponding backgate biases AVg 
measured from the Dirac point. 



localization length using the rough estimate [2£ 



~ Le exp 



e^/h 



(3) 



The elastic length Lg is derived from the Drude con- 
ductivity ao as Lf. = vpaD ■ /i/(2e^z)i?i/7m), assuming 
that iTD is the conductivity measured at a temperature of 
100 K. This yields £,d — 12 nm at the Dirac point, which 
is close to the variable range hopping estimate. The 
agreement is also quite reasonable at higher n (Fig. 6a). 
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Overall, this rough agreement gives us confidence in the 
order of magnitude of the localization length. 

The transition from weak localization to strong local- 
ization can be understood by contrasting the fundamen- 
tal transport length scales extracted from our experimen- 
tal data (Fig. 6a) . As long as the phase coherence length 
remains smaller than the estimated localization length 
5, the weak localization regime prevails. Whenever 
becomes comparable to ^, we observe that the conduc- 
tivity is close to e^//i, the value at which the transition to 
strong localization is expected Fig. 6a displays 
only when a > e^/h at zero magnetic field (and except 
at the Dirac point). In the opposite case (ct < e^//i), the 
comparison between experiment and weak-localization 
theory becomes worse, so that extracting a value for 
is meaningless. 

Graphene offers the possibility to tune the carrier den- 
sity n with Vg, which provides a practical knob to test 
the localization theory. From our measurements, we can 
construct a " phase diagram" of the transition from weak 
(WL) to strong localization (SL) as a function of n and 
temperature (Fig. 6b). To do this, we use the temper- 
ature dependence of u (a few such curves are shown in 
Fig. 3(b)) and define the transition as cr = /h. We em- 
phasize that the transition is expected to be gradual and 
to develop as a smooth crossover [s^. Figure 6b shows 
that the WL-SL transition is very sensitive to n at low 
carrier concentration. This is because of the strong vari- 
ation of ^ upon varying n (Fig. 6a). By contrast, the 
other fundamental transport lengths, Lg and L^, remain 
essentially constant. 

Eventually, one interesting outcome of our measure- 
ments is that they reveal the importance of Coulomb 
interaction between charge carriers. This can be seen in 
the magnetoresistance measurements at high magnetic 
field (Fig. 4a) where the conductivity at i? = 4 T is defi- 
nitely temperature dependent. The contribution of weak 
localization is reduced at high magnetic field and the re- 
maining correction to a is usually attributed to Coulomb 
interaction. The correction to the conductivity due to 
Coulomb interaction in a two-dimensional metal reads 
21, 



■ g2D ■ In 



(4) 



where g2D is a constant of the order of unity. Figure 7 
shows that experiment agrees well with the theory us- 
ing g2D — 0.8. Unlike weak localization, which origi- 
nates from a change in the diffusion constant 2l|, this 



Coulomb interaction effect is a correction to the density 
of states. These results indicate that the correction to 
the conductivity due to Coulomb interaction cannot be 
neglected in the transition between weak localization and 
strong localization in graphene, as it is treated in existing 
theoretical works [2|-l7|, 1321434 1 . 
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FIG. 7: Temperature dependence of the conductivity at B 
4 T (dots) and fit to equation (4) (solid line). 



CONCLUSION 

We have reported on the crossover from weak local- 
ization to strong localization in disordered graphene, the 
disorder being created with ozone. For this, we have car- 
ried out magneto-transport measurements as a function 
of gate voltage and temperature. We have shown that 
the transition between weak and strong localization oc- 
curs as the phase coherence length becomes comparable 
to the localization length and that the transition is very 
sensitive to the charge carrier concentration. In addition, 
we have demonstrated the importance of the resistivity 
correction due to disorder-induced electron-electron in- 
teraction. 

Previous works showed that disorder in graphene ob- 
tained b y hy drogenation or oxidation can open an energy 
gap 17h19{. By contrast, recent calculations reported 
that this is not the case for graphene exposed to ozone 
[3^ . In addition, even if such a gap was created, our mea- 
surements suggest that it would be very small. Indeed, 
we obtain an energy gap of 1 meV when we compare the 
temperature dependence of the resistivity at the Dirac 
point to a thermal activation behavior. This low value 
suggests that a gap, if it exists, would be relevant only to 
measurements near the Dirac point (in our case, it would 
be included in the interpretation of the measurements 
at = 36 V). Further experimental work is needed to 



demonstrate whether this energy gap exists or not, using 
e.g. scanning tunneling microscopy techniques. 
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